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The average transverse momentum transfer per unit path length to a fast parton scattering 
elastically in a perturbative quark-gluon plasma is related to the radiative energy loss of the parton. 
We first calculate the momentum transfer coefficient q in terms of a classical Langevin problem 
and then define it quantum-mechanically through scattering matrix element. After treating the 
well known case of a quark-gluon plasma in equilibrium we consider an off-equilibrium unstable 
plasma. As a specific example, we treat the two-stream plasma with unstable modes of longitudinal 
chromoelectric field. In the presence of the instabilities, q is shown to exponentially grow in time. 



I. INTRODUCTION 

When a highly energetic parton travels through dense 
QCD matter, it receives random kicks from elastic in- 
teractions with constituents of the plasma. The average 
transverse momentum transfer per unit path length is re- 
lated to the radiative energy loss of the parton [1| . The 
parameter describing the average amount of transverse 
momentum broadening per unit length is called q and is 
defined as 



q = d(Ap^)/dz, 



(1) 



when the fast parton flics along the direction z. The 
values of q extracted from experimental data on rela- 
tivistic heavy-ion collisions vary in a rather broad range 
0.5 — 15 GeV 2 /fm depending on the model of hard parti- 
cle propagation in strongly interacting matter produced 
in nuclear collisions [Hyj- 

The calculations of q for the case of perturbative quark- 
gluon plasma in equilibrium arc well understood (see 
[J, l5j for recent work). For such a plasma the value of q 
is predicted to lie at the lower end of the range of val- 
ues deduced from experiments 0]. However, it is not 
at all obvious that the quark-gluon plasma produced in 
a heavy ion collision is in equilibrium during the whole 
time interval of the propagation of the fast parton. It 
is rather expected that the momentum distribution of 
plasma partons is anisotropic at the early stage of rel- 
ativistic heavy-ion. And recently q has been computed 
[(| 0] for the case of homogenous quark-gluon plasma 
with anisotropic momentum distribution. However, such 
a plasma is unstable due to the chromomagnetic modes 
(for a review see Q). The fact that unstable systems are 
explicitly time dependent as unstable modes exponen- 
tially grow in time has not been taken into account in the 
previous analytical calculations 0, 0]- However, numer- 
ical simulations 0, [l(| indicate that q receives a sizable 



contribution from these unstable modes. Also, possible 
phcnomcnological consequences of the momentum broad- 
ening effect of an unstable plasma on hard partons have 
been proposed 0, [ll[ . 

In this paper we develop the approach suggested in 
[l2| (see also [l3[) to study the transverse momentum 
fluctuations in terms of classical Langevin problem. Wc 
calculate q by treating the parton as an energetic classi- 
cal particle with SU(3) color charge moving in the pres- 
ence of the fluctuating color fields. Then, q is expressed 
through the correlation function of chromodynamic fields 
computed in or equivalently, through by the un- 

ordered gluon propagator in Hard Thermal Loop (HTL) 
approximation. For the equilibrium plasma the Langevin 
approach recovers the known result obtained within the 
standard thermal field theory [l^|. We then apply the 
Langevin approach to the case of a parton traversing an 
off-equilibrium unstable plasma. For the sake of simplic- 
ity, we study the two-stream system and consider only 
the contribution from the unstable modes of longitudi- 
nal chromoelectric field, neglecting the contribution from 
transverse fields. Wc show that q grows exponentially in 
time due to the exponential growth of the unstable mode 
amplitude. 



II. CLASSICAL APPROACH 

We consider a classical parton which moves across a 
quark-gluon plasma. Its motion is described by the Wong 
equations [16j 



dr 



dp^{r) 
dr 



w m (t), 



9Q a (T)Fr(x(r))u v (r), 



(2) 



(3) 



2 



dQajr) 
dr 



= -gf abc p,(r)A^x(T))Q c (r), 



(4) 



where r, x m (t), u m (t) and p^{r) are, respectively, the 
parton's proper time, its trajectory, four-velocity and 
four-momentum; and denote the chromodynamic 
field strength tensor and four-potential, respectively, and 
Q a is the classical color charge of the parton. 

We look for a solution of the Wong equations in a spe- 
cific gauge and the gauge dependence of our results is 
discussed a posteriori. We assume that the potential 
vanishes along the parton's trajectory by imposing the 
axial gauge condition p^ (r) A^(x(t)) = 0. Then, Eq. Q 
simply states that Q a is constant as a function of r. 

We solve Eqs. ^ \5§ assuming that the momentum of 
the parton is so high that the changes in its direction 
of motion caused by the interactions with the medium 
are small compared with the momentum of the parton: 
I^pI/IpI <C 1. The change in the velocity vector v of 
the parton are then negligible, and we can consider the 
parton to move along a straight-line path with constant 
velocity. If the parton is massless or nearly massless (a 
light quark or a gluon) , its speed can be approximated by 
the speed of light. Although we assume further on that 
the fast parton is at least approximately massless, we 
note that the approximation of a straight line trajectory 
with constant velocity also applies to a heavy quark with 
mass M and momentum such that |p| = 7M|v| ^ |Ap|. 

In the limit of massless particle, the proper time be- 
comes a meaningless concept and we formally solve the 
equation of motion ^ to eliminate it. Then, we get 

l?(T)=p> t (0)+gQa f dr'F^{x{T'))u v {r'). (5) 
Jo 

In general, this is not an explicit solution of §3§ as 
the four- velocity u v (t) on the right-hand side of © 
depends on the four-momentum. However, the four- 
velocity, as explained above, is approximated by a con- 
stant: u v (t) w u v (0). Assuming that the direction of 
motion of the parton is in the positive x 3 -direction, we 
have explicitly u^ L = 7(1,0,0, v) with the Lorentz factor 
7 = (1 — w 2 )" 1 / 2 . Then (O becomes 

p»(t)=p»(0)+gQ a f dt'(Ff(x')-vFf(x')), (6) 
Jo 

where x' = x(t') and t = 7T. 

In the limit v — > 1, the form of © motivates the in- 
troduction of light-cone variables defined as 



P 



p°±p 3 



(7) 



In the new coordinates the four-vector components are 
understood as p^ = {p + , p~ , p 1 , p 2 ) and the product of 
two four-vectors is 



Assuming that the parton's initial position is x(0) = 0, 
its trajectory is given by x^(t) = s/2(t, 0, 0, 0), and 
Eq. ([6]) takes the form 



p"(z+) =p"(0) +gQ a I dy + Fr(y + ) 



(9) 



Here F£~(y + ) is a short-hand notation for F£~ (x l/ (y + )) 
with x v {y + ) = (y + , 0,0,0). If the parton moves in the 
negative cc 3 -direction, Eq. §§§ changes to 

p"(aT) =P fi (0)+gQ a [ X dy-F>?+(y-). (10) 



In the following we shall assume that the parton travels 
in the positive x 3 -direction. 

In the spirit of the Langevin approach, wc consider the 
ensemble average (p t *(x + )p L '(x + )) indicated by the angu- 
lar brackets. The ensemble average involves averaging 
over color charges. We recall that 



dQQ a Q b = C 2 5 a \ 



(11) 



where C 2 = 1/2 for particles (quarks) in fundamental 
representation of the SU(A C ) group and C 2 = N c for 
particles (gluons) in adjoint representation. Then, wc 
find 

(^(x+K(x+)) =p"(oy(o) 
,*+ 

+ g 2 C dytdy+ (Fr(yt)K-(vt)}, (12) 



(13) 



where the color factor C is given as 
C = 



jt- = jrr- (fundamental rep.) 



IvSri = (adjoint rep.). 



p IJ, q^=p T q +p <f -pV -p'q 



-,2„2 



(8) 



It is related to the eigenvalue Cr (R = F, A) of the 
quadratic Casimir operator as 

C R = G(N 2 C - 1). (14) 

Let us consider (p IJ, (x + )p fl (x + )} in more detail: 

{p^x+)p,{x + )) =p"(0)]v(0) 

+ g 2 C / dyf dyi g^K~ {vt)K~ ivt))- 
Jo 

(15) 

An interesting feature of (fT5|) is that although /1 runs 
over all four components 0,1,2,3 or +,—,1,2, only 
the transverse components (p = 1,2) contribute to 

g^( F a~(yi) F a~(yt))- This is easil y seen b y writing 

the expression out in light-cone coordinates and recall- 
ing that the antisymmetry of F" 13 requires F~~ = 0. 
The terms for [i = +, — thus vanish and consequently, 

g^-{y{)F:-{y 2 ) = 

- F^(yi)F^(y 2 ) - F 2 -( yi )F 2 -(y 2 ). (16) 
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The fact that only the transverse components (/x = 1,2) 
contribute implies that we can write 

(Ap|(.T+))^(Ap T (x+).A PT (x + )) 

= -g 2 C f dy+dy+ gilv {F>l-{yt)F»-{ y +)), 
Jo 

(17) 

where Apr(x + ) = px(x + ) — Pt(0). 

Equation (|17[) . which is the main result of this sec- 
tion, has been derived in a specific gauge and the right- 
hand side of the equation is, in general, gauge dependent. 
The gauge independence can be restored by inserting the 
Wilson line between the fields F£~(yf) and F^~(y2) in 
the correlation function (F£ ~ (yt)Fa~ (vt)) ■ However, 
the correlation functions, which are further used in our 
study, are computed in the Hard Loop approximation. 
These functions appear to be gauge independent even 
without the Wilson-line insertion, as discussed in detail 
in Sec. VIIIA of Q- Therefore, Eq. (fTT]) is gauge inde- 
pendent within the approximations we apply. 

We now assume that the quark-gluon plasma is trans- 
lationally invariant in space and time and, hence, the 
correlator (F£~ (y^)F^~ (y^)} depends only on the dif- 
ference y^ — y% . This assumption, which is relevant for 
equilibrium plasmas, will be not adopted for the two- 
stream system discussed in Sec. I VII as the growth of un- 
stable modes break the translational invariance in time. 
Making use of the translational invariance, we introduce 
the fluctuation spectrum 



/PC- 
a 1 a 



^(^(xjci, (is) 



which allows us to express Eq. p7[) as 



f x+ 

{Ap 2 T {x+)) = -g 2 C / dvtdyi 
Jo 



d^k 



e ik -^-^ gixv (F-»Fr)k. (19) 



Since the double integral over y^ and y£ tends to a delta 
function of k~ in the long-time limit 



dy+dy+ e lk (^"^) 
4sin 2 (fc-x+/2) x +_ 

we can define the transport coefficient 



2irx + S{k-), (20) 



q= lim -(Ap|(i)) = lim ^(Ap 2 T (x + )} 



1 

— i 



x+— too X" 1 



= -2V2ng 2 C J -0^ S(k~) g^Ff- F»~ ) k , 



(21) 



which was introduced in Eq. (fTJ as the transverse mo- 
mentum broadening per unit path length. 

From here on we can proceed in two different ways. 
We can either use the fluctuation spectra of chromody- 
namic fields derived in [l4|, or we can express (F£~ F% ) k 
through the unordered gluon propagator in the HTL ap- 
proximation. In the two following sections we pursue 
both approaches. The first one is shorter, while the sec- 
ond one allows us to relate our results to those obtained 
in0. 



III. FLUCTUATION SPECTRA OF 
CHROMODYNAMIC FIELDS 

The fluctuation spectra of the equilibrium plasma, as 
computed in [bH ]. are 



{E l a El) k = 28 ab T^ 



k l ¥ 3e L (w,k) 
k 2 | W 2 £ L (u;,k)| 2 
k l W\ 3e T (cj,k) 



(B % a B 3 b ) k = 2<5 ab Tuk 2 (5^' 
3e T (w,k) 



k 2 J |w 2 e T (w,k)-k 2 | 2 



k 2 



|w 2 £ T (w,k)-k 2 | 2 ' 



{B\El) k = {ElBl) k 



25 ab Tuj 2 e lk] k k 
3£t(w, k) 



(22) 



|w 2 e T (w,k) -k 2 | 2 ' 

where k = (w,k), e 1 ^ is the antisymmetric tensor, T 
is the system's temperature and £z,,t(w, k) are choro- 
modielectric functions which for the equilibrium plasma 
of massless particles are well-known to be (l5| 



^e L (cu,k) 
Sy£i(w,k) 
^£ T (w,k) 



= 1 



k 2 





LU + 


|k| 




1 - -fLln 








. 2|k| 


to — 


|k| 





2k 2 



(23) 



2w k| 



-In 



9£ T (w,k) 
where ttld is the Debye mass 
„a 9 2 T 2 



^6(k 2 -. 2 ) W - (k2 - w2) 
4 



^|kl : 



' D 



G 



{N f + 2N C ). 



(24) 



(25) 



Since 



9,AFrK-)k = {E x a E:) k + {EyEt) k 
(E*B») k + (E»BZ) k 
{BlEl) k + (BlEl) k 
+ (B»B*) fc + {BlB») k , 



(26) 
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Equation (|2l"j) yields 



2g 2 C B T 



d 3 k 
$e L (k z ,k) 



k z k 2 

k\k\ 3e T (fc z 



|e L (fc z ,k)|2 |fc 2 e T (fc 2 ,k)-k 2 | 2 



(27) 



The classical formula (f27|) holds for the inverse wave 
vectors of the fields which are much longer than the 
de Broglie wavelength of plasma particles. It requires 
|k| <C T where T is the plasma temperature. For larger 
wave vectors a quantum approach is needed which is dis- 
cussed in Sec. [V] 

Keeping in mind that the equations 

£i(w,k)=0, Lu 2 e T {Lo,k) - k 2 = 0, 

determine the longitudinal and transverse collective 
modes in the isotropic plasma, one sees that the trans- 
verse momentum broadening, as given by the formula 
(|27| . is determined by the interaction of the fast parton 
with plasma collective excitations. 



IV. GLUON PROPAGATOR 

As already mentioned, (F£ ~F^~) k , which enters 
Eq. (f2"Tj) . can be expressed through the unordered gluon 
propagator in the HTL approximation. 

Neglecting all non-Abelian contributions, the Fourier 
transform of the field strength tensor is expressed 
through the four-potential as 



Fr(k) = -ik»A» a (k) + ik»AZ(k), 
allowing us to rewrite (j2~lj) as 



(28) 



where, as already mentioned, the color indices are sup- 
pressed. We also introduce the spectral function 

P^ix-y) = ([A^x),^)]) 

= iD> I/ (x-y)-iD< I/ (x-y). (32) 

The equilibrium Green functions D^ v (x), D< v (x) obey 
the Kubo-Martin-Schwingcr condition 



or 



D> (i,x) =D<(t + 0,x), 



D>(k) = e?» D<(k), 



(33) 



(34) 



where (3 = 1/T. Combining the definition (f3"2"| and the 
KMS condition lIMl). we have 



(35) 
(36) 



iD> u (k) = (n(w) + 1) /v(fc), 



where n(co) = (e 13 ^ — 1) _1 denotes the Bose distribution. 

The spectral function in the HTL approximation can 
be expressed as [l5| 



P^u{k) 



2 9eL ( W> k) pL 



fc 2 |e L (w,k)| 2 * 
2w 2 3e T (w,k) 



+ 



|w 2 e T (w,k)-k 2 | 2 " 
with the projectors defined as 



P^(k) = s ij - 



k l w 



pr(k) 



k 2 ' 

k^k v 



P^(k) = 0, 



P^{k). 



(37) 

(38) 
(39) 



q = V2irg 2 C J ^ 5(k~) fc 2 (A:A~) k . (29) 



Here we made use of the relation d°±d 3 = \^2d T and fur- 
ther used the fact that only the fx = 1, 2 components con- 
tribute in the field strength correlator. We also dropped 
all terms proportional to k~ in view of the delta function. 

Because y^ and y^ in Eq. (|17|) run independently from 
to x + , the correlator (A~A~) k is the Fourier transform 
of the time-unordered gluon Green function. Since the 
plasma of interest is on average color neutral, the gluon 
propagators such as (A£(x)A%(y)) are proportional to 
5 ab . We drop for now the color indices but when the 
gluon propagator (A~A~)k is substituted in Eq. (|2"9")l . 
the coefficient 8 aa — TV 2 — 1 will included. 

For translationally invariant systems, one defines the 
time-unordered (>,<) Green functions 



iD> v (x-y) = (Ap(x)A v (y)) t 
iD< v {x~y) = (A^A^x)), 



(30) 
(31) 



We identify the correlator (A a A b from Eq. (j2"Tj) with 
the time-unordered Green function D"(k). Substitut- 
ing D>~(k) given by Eqs. ([SSI [S3) into Eq. (gSJ, we ob- 
tain 



q = aV2tt 9 2 C r J S(k-) n(w) k% (40) 

%s L (oj,k) 



k 2 \e L (cu,k)\ 2 L 
cj 2 £t(w, k) — k 2 



P£-(k) 



p--(k) 



12 1 T 



where Cr is the eigenvalue of the quadratic Casimir op- 
erator given by Eq. ([33]). 

With the help of formulas ([351 ESI) we find 



P t (k) - 77i3' 



2k 2 : 
k~k~ 



k 2 



2k 2 ' 



(41) 
(42) 
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and thus, we finally obtain 



2g 2 C R 



d 3 k 



(2tt) 3 
%e L (k z ,k) 



k 2 

»(*.) § 



\a L (k z ,k)\* \k 2 e T (k z ,k)-h 2 \ 2 



(43) 



where k z = k 3 . In the classical field limit, where n(u>) ~ 
T/w, we obtain again the expression (f2T)) . 

Since the HTL approximation holds for |k| <C T, the 
integration in Eq. (|4"3"| . as well as in the fully classical 
formula (j2~T)) . should be cut off at |k max | -c T. For larger 
wave vectors we need a quantum approach, which is dis- 
cussed in the next section. 



V. QUANTUM APPROACH 

Here we sketch the derivation of q following the strat- 
egy outlined in [l5[ for the electromagnetic energy loss 
of a fast muon. We present how to obtain two contribu- 
tions to q coming from soft (|q| <C T) and hard (|q| ~ T) 
momentum transfers. We discuss only the main steps of 
the derivation as the procedure is mostly described in the 
literature 0, 0| ■ This section is included for completeness 
of our study. 

Within a quantum approach q is defined as 



j2 2 

d q T q T 



d 2 T 

d 2 q T 



(44) 



where T is the rate (probability per unit time) for elas- 
tic collisions of parton with plasma particles. For a fast 
parton with four- momentum p\ which scatters on plasma 
constituents of species i with four- momenta p 2 , the rate 
corresponding to the binary process p\ + P2 — > p[ + p' 2 is 



given as 



r = — 



d 3 p[ d 3 p' 2 d 3 p 2 



(45) 



2E' 1 (2n) 3 2E 2 (2n) 3 E 2 (2n) 3 
x(27r) i 5^(p 1 +p 2 -p' 1 -p' 2 ) 

x Yj ( p2 ) t 1 ± /i(Pa)] \Mi(pi , P2;p[,P 2 )\ 2 , 

i 

where /i(p) is the distribution function of particles of 
species i and A4i(pi,p 2 ;p' 1 ,p 2 ) is the scattering matrix 

I 1 2 

element. The definition of \M-i\ involves summation 
over final spins and colors and an average over the initial 
spin and color of the fast parton. 

The matrix element of parton-parton scattering is as- 
sumed to be dominated by the onc-gluon exchange pro- 
cess. Then, the matrix element squared can be expressed 
as the self energy of the fast quark. When the fast parton 
is a massless quark, one finds Q 



g=^Tr[/5£>(p)] 



(46) 



where S^(p) is the fast quark self-energy 'weighted' with 
q\. It is given as 

K(P) = 9 2 C F J 2£ ^ )3 ^W(g), (47) 

where p' = p + q and D< u (q) is the gluon propagator 
(stripped of the color factor 5 ab ) which includes the one- 
loop correction. Changing the orientation of the mo- 
mentum q, the propagator D< v (q) ([55]) should be re- 
placed by D> v (q) |36|) but we get the same result as 
n{-qo) + 1 = -n{qo). 

In the soft domain, where |q| <T, the one- loop gluon 
propagator D<„(k) is obtained in the HTL approxima- 
tion. Substituting D< v (k) provided by Eqs. ([331 IB"?)) 
into Eq. (|46l) . we find the differential rate which enters 
Eq. (01) aTQj 



d 2 T 



soft 



d 2 q T 



= 2g 2 C h 



dq z n{q z ) 



(2tt)S 



^se L (q z ,q) q\q\ 5s T (q z , q) 
Mfe,q)| 2 + |^T(9-.k)-q 2 | a 



(48) 



where the upper cut-off A obeys gT < A < T. We have 
here replaced Cp by Cr as a similar analysis can be per- 
formed for fast gluons. The expression for q given by 
Eq. ((33]) with (gSJ) agrees, as expected, with Eq. 1(15]). 
By explicit evaluation one shows [J] that the soft colli- 
sion rate ()48|) can be approximated as 



dT 



soft 



Cr 



Tm 2 D 



d 2 qi 



{2w) 2 q 2 {q 2 +m 2 D )' 



(49) 



The hard (|q| ~ T) contribution to q can be found di- 
rectly from Eqs. (|44|45|) using the unscreened gluon prop- 
agator. The approximate result is [H 



dT 



hard 



d 2 qi 



g 



Cr p 
(2ir) 2 q\ 



(50) 



where p is the properly normalized parton density in the 
plasma. 

The complete result for q is obtained by summing up 
the hard and soft contributions. An arbitrary parame- 
ter A, which separates the soft momentum transfers from 
the hard one, is the ultraviolet cut-off for the soft con- 
tribution and the infrared cut-off for the hard contribu- 
tion. A is eliminated from the final answer, as the soft 
contribution depends logarithmically on A and the hard 
contribution depends logarithmically on A -1 . Since the 
hard contribution also depends logarithmically on the ul- 
traviolet bound provided by two-body kinematics of the 
scattering process, we can simply extrapolate the soft 
contribution to the hard domain with no significant er- 
ror Q . In conclusion, the integral in the soft contribution 
to q, as given by Eqs. (|44I49|) . can be extended to the 
maximal momentum transfer g™ ax or fc™ ax roughly equal 
\J ET where E = p° is the energy of the fast parton. 
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VI. TWO-STREAM UNSTABLE PLASMA 

Our aim here is to calculate the momentum broaden- 
ing of a fast parton in unstable anisotropic plasmas. As 
mentioned in the Introduction, this problem was studied 
previously @, 0] in a way similar to that presented in 
Sec. |Vl Namely, the hard thermal loop gluon propagator 
of an equilibrium plasma was replaced by the hard loop 
propagator of an anisotropic plasma [181 ]. Such an ap- 
proach, however, does not take into account that unstable 
systems are intrinsically time dependent because the un- 
stable modes grow exponentially in time. We show here 
how these modes influence the momentum broadening of 
a fast parton. For the sake of analytical tractability, we 
do not consider the anisotropic plasma with a momentum 
distribution likely to be found in the early stage of rela- 



tivistic heavy-ion collisions, as the correlation functions 
of chromodynamic fields are difficult to obtain for such a 
plasma. Instead, we consider a fast parton in the more 
tractable, but also unstable two-stream plasma. With 
this simplified example we hope to elucidate the general 
features of the problem of momentum broadening in un- 
stable plasmas. 

The distribution function of the two-stream plasma is 
chosen to be of the form 

/(p) = (27r) 3 n [S^ (p - p ) + 6& ( p + Po) ] , (51) 

where n is the effective parton density in a single stream. 
The distribution function (|5ip should be treated as an 
idealization of the two-peak distribution where the par- 
ticles have momenta close to po or — po- 



J 



The two-stream plasma is unstable with respect to both longitudinal and transverse modes. Because the two- 
stream plasma does not have a direct phenomenological application, and we are only interested in showing the effect 
of unstable modes on q in principle, we consider for simplicity that the parton interacts only with the longitudinal 
electric fields. It is worth noting that there exist also unstable transverse modes, which will have a similar effect on 
q. The correlation function of the longitudinal field components was derived in 14|. We are mostly interested in 



the exponentially growing modes but the correlation function which represents only such modes becomes meaningless 
when the growth rate tends to zero. Therefore, we denote the correlation function with the subscript label "cxp" to 
indicate that it includes both, the exponentially growing and the exponentially decaying modes. The expression for 
the correlation function reads UM '■ 



(^(t liri )^(i 2 ,r 2 )) 



5 ab n 



d 3 k e ik( - ri - r2) 



k l W 



U+jk-u) 2 )' 



2 " "J (2^) 3 k 4 (lU 2 + -LU 2 _) 2 

7k + (k • u) 2 ) cosh (7,^ + i 2 )) + ( 7 £ - (k • u) 2 ) cosh ( 7k (ti - i 2 )) 



(52) 



where u = po/-B Po is the stream velocity and ±w±(k) are the four roots of the dispersion equation ^(t^k) = 0: 



4(k) = V2 



k 2 (k • u) 2 + ^ 2 (k 2 - (k • u) 2 ) ± n J(k 2 - (k • u) 2 ) (4k 2 (k • u) 2 + M 2 (k 2 - (k ■ 



(53) 



with fi 2 = g 2 n/2E po . 7k is the instability growth rate 
(0 < 7k € R) defined as w~(k) = i7 k . The integration 
over k in Eq. (|52p is performed over the domain of un- 
stable longitudinal modes defined by the requirement 



k-u^0 and k 2 (k-u) 2 < 2n 2 Ck 2 - (k ■ 



(54) 



The correlation function (|52| is obviously invariant with 
respect to space translations - it depends on the differ- 
ence (ri — r 2 ) only. The initial plasma state is on average 
homogeneous and it remains so over the course of time. 
The time dependence of the correlation function (|52|) is 



very different from the space dependence, and it is not in- 
variant under time translations. The electric fields grow 
exponentially and so does the correlation function, both 
in the variables (ii + < 2 ) and (t\ — 1 2 ). The fluctuation 
spectrum also evolves in time as the growth rate of unsta- 
ble modes is wave vector dependent. After a sufficiently 
long time the fluctuation spectrum is dominated by the 
fastest growing modes. 



Substituting the correlation function ([52"]) into Eq. (fT5|) 
and reverting to integration over time t instead of the 
light-cone variable y + , we find 
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'o Jo 



(ApUt)) = g 2 C dt> / dt" (EZ{t\T{t')EZ{t>\r{t")) + {El{t',T{t')El{t",r{t")) 



( 7 2 + (k-u) 2 y 



(<£ + (k • U) 2 ) (| e (*--HTk)* - 1| 2 + | e (*--7k)t _ 1)2 ) + 4 ^2 _ (k . u) 2) ; 



7 k 



7£. 



Taking into account only the fastest growing contribution to Eq. (|55p. we obtain 



(55) 



„ rf(Apf^)) <? 4 _ /• d 3 fc 

9 ~ — * y CRn y 



>27k* 



4( 7 2 + (k.u) 2 



k^ 2 - W 2 ) 2 7k (fc 2 + 7 2 )' 



(56) 



This expression diverges when 7 k — > 0. It happens be- 
cause only the unstable modes are taken into account. 
The formula given in Eq. (|55[) remains finite in the limit 
7 k — > 0. The momentum broadening (|56|l grows exponen- 
tially in time as the spontaneously growing fields exert an 
exponentially growing influence on the propagatingpar- 
ton. This effect is missing in the previous results [g, 0] 
obtained for the unstable anisotropic plasma, which was 
treated as a stationary medium. 

The calculation of the momentum broadening carried 
out here does not assume that the medium is close to 
equilibrium or that the growth rate of the instabilities is 
small enough to invoke the static approximation. Even 
though we have restricted our calculation to the simple 
and analytically tractable case of the two stream system, 
the underlying formalism may be applied to calculate the 
momentum broadening in an arbitrary medium as long 
as the field strength correlators are calculable. 

The differences between the results obtained in this 
manuscript and those in previous articles using a steady 
state approach are two-fold: An unstable system is not 
time-translation invariant and thus a calculation based 
on time translation invariancc is in general invalid. When 
the propagator depends independently on t and t' , rather 
than on the difference t — t' , one obtains a function of 
two frequency variables when performing the (one-sided) 
Fourier transformations with respect to the two time vari- 
ables. 

In assuming that the gluon occupation number n(u)) 
is a function of only a single frequency, authors of 
Rcfs. 0, 0] implicitly assume that the rate of instability 
growth is slow enough so that one may impose approxi- 



mate time translation invariance over certain time peri- 
ods. Note, however, that even if we ignore the {t\ + t 2 )- 
dependence in the first term of the correlator (j52")) . the 
Fourier integral with respect to t\ — t<i of the second term 
does not exist because of the exponential divergence. If 
we were to ignore this problem for the moment, we would 
obtain results which could be reconciled with those of 
Refs. only if we would make the additional assump- 
tion that the system would remain very close to thermal 
equilibrium. In the following, we illustrate this point in 
more detail by contrasting our calculation with the re- 
sults obtained in the steady-state approach for the two- 
stream system. It is useful to start from the expression 
(|43|) for the parameter q. 

In order to decompose the expression for q into a con- 
volution of an occupation number and a spectral density 
we now assume approximate time translation invariancc. 
We again restrict our consideration to the longitudinal 
response of the plasma. The calculation is most easily 
carried out in the temporal axial gauge, in which the 
spectral function of longitudinal electric fields is given 
by 



2 3e L (fc) k l ti 



\e L {k)\ 2 k 2 



(57) 



A straightforward calculation shows that this expression 
continues to hold when the system is anisotropic, but 
the anisotropy is determined by a single vector, here the 
flow velocity u. The dielectric function £x,(o;,k) for the 
two-stream system is given in Ref. [l4j : 



£L(w,k) = 1 - 



(k-u) 5 



k 2 



(w-k-u) 2 (w + k-u) 2 



(^ 2 -^ + (k) 2 )(co 2 -cMk) 2 ) 
(w 2 - (k-u) 2 ) 2 



(58) 



where w±(k) is given by (|53[) . We note that the expression for 3£l(w, k) vanishes for real u> and k for the two-stream 
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system. The expression SJex(o;, ]s)/\ei,{u>, k)| 2 in (|57|l should therefore be replaced by 



1 



e L (w,k) + isgn(w)0 H 



7rsgn(w) <5(ei(w,k)) 



(59) 



The expression analogous to (|55|) derived from the steady-state formula then reads: 



d 4 ifc fc 2 



e L (£; z ,k) + isgn(fc z )0 H 



<i 3 /fc A: 2 
■Kg 2 C R I /r> _ vc , rt(fc 2 ,k) ^ sgn(fc 2 ) (5(£ £ (fc z ,k)) 



(2tt) s 



(60) 



Here n(aj, k) is the average occupation number of the soft gluon modes. For the case of a medium in thermal 
equilibrium, n{uj, k) is given by the Bose distribution and for soft frequencies may be approximated by T/lo. For 
the case of a medium that is barely anisotropic, i. e., where the corrections to the Bose distribution are small, one 
may persist with this form for the occupation number for short periods of time. This is the approximation made in 
Refs. 0, 0]. Such an approximation will, no doubt, break down as the medium continues to depart from isotropic 
equilibrium or as the unstable modes grow with time. Besides the insistence on the static approximation, the choice of 
the equilibrium Bose distribution without specifying the short range of time over which this is applicable, constitutes 
an additional source of error in these earlier works. ^ 



In order to obtain the correct form of the occupation 
number for any time, in the limit that one may still invoke 
the static approximation, one needs to recalculate the two 
unordered space-time propagators, 



D>(t) = Y^(n\pA(t)A(0)\i 

It 

D<(t) = £<n|M(0M(f)|» 



(61) 



for the specific density operator characteristic to the un- 
stable system being considered. In the equation above, 
A(t) is the vector potential of the glue field suppress- 
ing Lorentz and color indices, as well as space variables. 
In this case one may obtain the spectral density (|32p in 
momentum space and thus calculate the time dependent 
average occupation number as 

D<(w) 



(62) 



D>(lu) - D<{uj) 

Only in the case that the density operator is given by 
the thermal Gibbs expression p = cxp[— f3H], will there 
be a simple Kubo-Martin-Schwinger (KMS) relation (f34|) 
resulting in the Bose distribution for the occupation num- 
ber. For any other form of the density operator, the oc- 
cupation number may be considerably different from the 
Bose distribution, even under the assumption of approx- 
imate time translation invariance. 

As a result, we find that the use of a static approxi- 
mation and a Bose distribution for the occupation num- 
ber of the unstable modes restricts the applicability of 
the calculations of Ref. @, 0] to very short times in a 
plasma where the growth rates of the instabilities are 
rather small. It should come as no surprise that in such 
a regime, the effect of such modes on the momentum 
broadening of a hard parton is minimal. In an arbitrary 
weakly coupled quark-gluon plasma with a large momen- 
tum anisotropy, the growth rate of the unstable modes 



may be large and, as a result, the calculations in the cur- 
rent article remain the only consistent approach to cal- 
culating the broadening of a hard parton traversing such 
a system. Furthermore, wc note that the badly singular 
contributions encountered by the authors of Ref. [g, 0| 
in their calculations presumably arise from the improper 
formulation of the problem as explained above. 



VII. SUMMARY AND OUTLOOK 

We have studied the effect of color field fluctuations 
in a QCD medium on the propagation of a hard parton. 
In the first part of this work, we described a classical 
Langevin approach to compute the transport coefficient 
q which denotes the amount of transverse momentum 
broadening per unit path length of an energetic parton. 
In this Langevin approach the coefficient is expressed 
through the correlation function of chromodynamic fields 
along the light-cone or, equivalcntly, through the un- 
ordered gluon propagator. In the case of a perturbative 
equilibrium plasma when the field correlation functions 
or the gluon propagators are described in the HTL ap- 
proximation we reproduced the known result 0, @] • 

In the second part of this study we considered a hard 
parton in an off-equilibrium, unstable plasma. Our ex- 
plicit calculations performed for the two-stream plasma 
showed that unstable modes with exponentially growing 
amplitude cause an exponential growth of the transverse 
momentum broadening. Our result demonstrates that 
the quasi-static approximation used in previous calcula- 
tions of the transport coefficient q is not applicable. 

It should be clearly stated, however, that the temporal 
evolution of anisotropic plasma is not the only feature 
which needs to be taken into account when transport 
properties of the unstable plasma are studied. When 
the plasma is populated with strong fields, new inter- 
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action mechanisms appear. For example, it has been re- 
cently argued that in a spontaneously chromomagnetized 
plasma the viscosity acquires an anomalous contribution 
[19( while synchrotron radiation contributes to the par- 
ton energy loss 0, Hl| . 

The formalism for the calculation of q used in this work 
is applicable to numerical simulations. It would be inter- 
esting to perform such a calculation for the case of a 
realistic, three-dimensional quark-gluon plasma off equi- 
librium. We note that the calculations of the field corre- 
lators are very complex even in the linear response regime 
[HI . When the growth of a continuum of unstable modes 
results in strong nonlinear interactions among them, as 
demonstrated in the case of the anisotropic plasma by 
numerical simulations p2l . [23[, an analytical treatment 
of the field correlators seems to be impossible. In view 
of the great phenomcnological importance of the trans- 
port coefficient q in relativistic heavy ion physics, a nu- 
merical evaluation of the field correlator in an unstable, 
anisotropic plasma would be very desirable. According 
to the exploratory simulations 0, U3] , q can be signifi- 
cantly larger in the unstable plasma than in an equilib- 
rium plasma. 

Another promising direction for future study is the ex- 
tension of our formalism to momentum broadening in a 
strongly coupled plasma. Due to the asymptotic freedom 



of quantum chromodynamics, a highly virtual, energetic 
parton can be weakly coupled to the quark-gluon plasma, 
even if the dynamics of the plasma itself is governed by 
strong coupling. Such an investigation would, for exam- 
ple, pursue the question of potentially large logarithmic 
corrections to q arising from large ratios of scales, such as 
Q 2 /T 2 , where Q 2 denotes the virtuality of the hard par- 
ton. Instead of relying on an expansion in orders of the 
coupling constant g, an approach based on factorization 
into nonperturbative matrix elements of the medium and 
perturbative dynamics of the parton itself would lend it- 
self to a renormalization group treatment [24| ■ Although 
q would not be calculable ab initio in such an approach, 
its scale evolution with the energy or virtuality of the fast 
parton could be predicted. 
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